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Abstract 


Owing to the immobility of traditional reactors and spallation neutron sources, the demand 
for compact thermal neutron radiography (CTNR) based on accelerator neutron sources has 
rapidly increased in industrial applications. Recently, thermal neutron radiography experiments 
based on a D-T neutron generator performed by Hefei Institutes of Physical Science indicated a 
significant resolution deviation between the experimental results and the values calculated using 
the traditional resolution model. The experimental result was up to 23% lower than the calculated 
result, which hinders the achievement of the design goal of a compact neutron radiography system. 
A GEANT4 Monte Carlo code was developed to simulate the CTNR process, aiming to identify 
the key factors leading to resolution deviation. The effects of a low collimation ratio and 
high-energy neutrons were analyzed based on the neutron beam environment of the CTNR system. 
The results showed that the deviation was primarily caused by geometric distortion at low 
collimation ratios and radiation noise induced by highenergy neutrons. Additionally, the 
theoretical model was modified by considering the imaging position and radiation noise factors. 
The modified theoretical model was in good agreement with the experimental results, and the 
maximum deviation was reduced to 4.22%. This can be useful for the high-precision design of 
CTNR systems. 
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1 Introduction 


Neutron radiography has been established as an important non-destructive inspection method 
and quantitative measurement tool 1. Because thermal neutrons can provide a strong contrast for 
some elements close to one another in the periodic table and can even distinguish isotopes of the 
same element, neutron radiography can provide complementary X-ray and gamma-ray 
radiography "!. It is widely used in the aerospace industry, national defense, materials energy, 
biological archaeology, and other fields. Currently, thermal neutron radiography (TNR) is an 
important field in the development of neutron radiography PI. However, the traditional TNR is 
mainly based on an expensive, bulky, and immovable reactor or spallation neutron source, which 
greatly limits its application. Compact thermal neutron radiography (CTNR) can mitigate these 
challenges 151, 

In contrast to conventional TNR, CTNR commonly uses a compact accelerator as the neutron 


generator. Accelerator-driven neutron sources (such as D-D or D-T neutron generators) require 
complex moderator and collimator configurations to create neutron beams with suitable 
radiographic properties '!. Its intensity is 2-3 orders of magnitude lower than that of the reactor or 
spallation neutron source. Consequently, the quality of the available neutron beam after 
moderation and collimation is significantly lower than that of conventional neutron radiography, 
which mainly manifests as a low collimation ratio (L/D) (Table 1) and low thermal neutron 
content (Fig.1) 57]. Both of these factors have an adverse effect on imaging resolution, but are less 
considered for conventional TNR systems. Based on this background, this paper presents an 
in-depth study of the resolution of a CTNR system. 
Table 1 Parameters of TNR facility 


Type Neutron source Imaging L/D Imaging area Spatial resolution 
neutron flux (cm?) (um) 
(n/s/cm?) 
Reactor IBR-2 [8] ~5.52x10° 200-2000 20x20 300 
Reactor WWE-K "1 ~7.1x10° 75-1400 20x20 232.2 
Reactor PULSTAR [1 1.8x10°-7x 10° 100-150 2.7*2.7 84 
Spallation RANS H] ~104 21-781 17x17 40 
Accelerator BARC !?] ~4x10 20 20x20 ~1000 
Accelerator PKUNIFTY [1] 2.35104 50-200 21x21 330 
Accelerator CAEP l4 ~1x10+ ~25 20x20 500 
Accelerator INEST ~1x104 ~15 20x20 600 
_~ 10!° 10! 
a — Spectrum from light water reacto 
E Spectrum from moderated prg 
pe accelerator driven neutron source Ta 
= 10" 10° “g 
t : 
5 Š 
S 10'4 108 E 
5 Z 
© 
É 
a 108 daanin susud al suil sti 1 7 
108 107 10° 10° 104 10° 107 10' 10° 10! 


Energy (MeV) 


Fig.1 Comparison of neutron spectrum between the reactor and moderator 
2 Shortcoming of the traditional theoretical resolution model 


2.1 The traditional resolution model of neutron radiography 


The modular transfer function (MTF) is a classical resolution analysis method based on the 
frequency domain derived from the Fourier transform of a point spread function (PSF). Compared 
with the PSF and other spatial domain methods, the MTF method has higher accuracy in 
resolution calculation. The MTF is defined as the ratio of the output modulation to the input 
modulation. The MTF of an image can be calculated using the following equation !!>!: 


MTF(f) me M out T Crna! nn) | Cnatan) 
Ma, Onnaxtmin)/ Omas +Inin) 


(1) 


where I,a is the maximum grey value at the target area of the output image, T is the 


min 


minimum grey value of the output image, Imay is the maximum grey value of the target area of 


the input image, and J,,;, is the minimum grey value of the input image. According to the 
Rayleigh criterion, the spatial resolution of an image is defined as the corresponding resolution of 
the 10% MTF "'°l, Furthermore, the MTF method can be used to analyze the effects of system 
elements on the spatial resolution. 

The model of resolution is essential for designing a TNR system to ensure that realistic 
resolution goals are established and achieved. The traditional theoretical resolution model focuses 
on the TNR system structure, which comprises a collimator, converter screen, and imaging system 
USI. Assuming that the effects of these components on the imaging resolution are independent of 
each other, the total MTF of a digital thermal neutron imaging system can be analytically 
produced through a Fourier transformation !!71: 


MTF (u )=sine (su exp (-m (6u 
Y)y-sine Gis) y 


where u is the spatial frequency of imaging, d is the distance of the converter screen from the 
sample, D is the diameter of the neutron aperture, L is the distance from the aperture to the sample, 
ô is the optical diffusion response of the converter screen, As is the sampling at the sensor, Mccp is 
the scintillator-to-CCD magnification, and As / M cep is image sampling on the 
converter screen. The MTF describes the magnitude of the system’s frequency response. The 
theoretical resolution model is helpful for depicting and quantifying the system resolution. 


2.2 The deviation of the theoretical resolution model applied to CTNR 


Based on the traditional theoretical resolution model, this study analyzed the factors affecting 
the imaging resolution of a TNR system driven by the typical accelerator neutron source listed in 
Table 1. Fig. 2 shows the relationship between the resolution and L/D with different optical 
diffusion responses and image samplings on the converter screen. When the L/D ratio is less than 
12, the curves for each parameter coincide. This result indicates that the influence of the converter 
screen and image sampling is small in the resolution analysis at a low L/D ratio and that the L/D 
ratio becomes the most important factor affecting the resolution. 
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Fig.2 Relationship between resolution and L/D: (a) optical diffusion response and (b)image sampling. 

To verify the applicability of the theoretical model, a standard line-pair (Fig. 3b) sample 
neutron radiography experiment was performed based on the TNR terminal with a D-T fusion 
neutron source (Fig. 3a) built by the Hefei Institute of Physical Science (HFIPS), Chinese 
Academy of Sciences. The iron sample had dimensions of 50 mm x 100 mm x 10 mm and seven 
line pairs of 0.5, 0.62, 0.83, 1.0, 1.67, and 2.0 lp/mm separately. A schematic of the system is 
shown in Fig. 3c. Table 2 lists the characteristics of the neutron beam emitted by the moderated 
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collimator of the neutron radiography terminal and the information of the imaging system. 


Accelerator 
D-T neutron Be 


Converter 


Screen Mirror 
Sample ae 


Collimator 


Accelerator 
D-T neutron ; 
source — -- aMis- 


Light 
tight 
Box 


Moderator 


Output 


Fig.3 Experimental setup for the HFIPS: (a) picture of the terminal, (b) picture of the line-pair sample, and (c) 
schematic of the setup. 


Table 2 Parameters of compact thermal neutron radiography facility in INEST 


Parameter Value 
Thermal neutron flux 1.5x10° to 2x104 n/cm?/s 
L/D 7.5 to 15 
Divergence ~ 85° 
Field of View (FOV) 20 x 20 cm? 
Converter screen SLiF/ZnS (Ag) 1:2 
Pixel number 1024 x 1024 pixel 
Pixel size 13 um x 13 um 


Neutron radiography experiments were conducted using L/D values of 6.3, 7.4, 8.8, 10, and 
12, respectively. The imaging exposure time was 300 s. Fig. 4 shows the radiography results of 
line pairs with different L/D ratios. The center of the imaging field of view (FOV) coincided with 
the center of the neutron beam. Fig. 4a shows the neutron radiograph with L/D=7.4. The red box 
in the figure represents the sampling area of the line-pair image, and the corresponding gray curve 
was obtained. The gray curves of the images with other L/D ratios are shown in Fig. 4b. The 
corresponding MTFf) function value of each line-pair image can be obtained by substituting the 
gray value into Eq. 1. Then, the image resolution was obtained according to MTF=0.1 by using the 
interpolation method, as shown in Fig. 5a. The theoretical resolution results under the same 
conditions were calculated using Eq. 2. A comparison between the experimental results and the 
theoretical resolution curve is shown in Fig. 5b. The results show that the image resolution 
improves with an increase in L/D in both the experimental data and theoretical calculations. 
However, all the experimental points were higher than the theoretical calculation curve, and the 
maximum deviation reached 23.1%. There was a significant deviation between the theoretical 
model and the experimental results, and the reason for this deserves further study. 


©) L/D=6.3 Tyyax=105.1 Iy,,=40.3 L/D=8.8 Tyggx=143-1 Tyjg=80 
L/D=10 T. =165.1 In„=22.3 L/D=12 OT =82.5 Imi =63 


Fig.4 Results of the standard line-pair sample experiment: (a) result with L/D=7.4 and corresponding sampling area; 


(b) results with other L/D ratios. 
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Fig.5 MTF Calculation results: (a) under different L/D conditions, and (b) comparison between experimental results 


and theoretical curves. 


3 Analysis of main parameters affecting spatial resolution and 


resolution deviation 


In this section, the effects of a low L/D ratio and a high proportion of nonthermal neutrons on 
the imaging resolution are discussed to elucidate the reasons for the deviation mentioned above 
based on the Monte Carlo simulation numerical analysis method, which is performed on the 
GEANT4 framework !!3l, The G4TENDL database [!°l, which has been proven to have high 
computational accuracy in previous studies ?°!, was used in the Monte Carlo calculations. 


3.1 Effect of low L/D ratio on CTNR system 


3.1.1 Geometric distortion 


In this study, a GEANT4 Monte Carlo calculation model was built based on a theoretical 
resolution model. As shown in Fig. 6, the model consisted of three modules: a surface neutron 
source with a diameter of 9 cm, a converter screen with an area of 20 x 20 cm’, and a 
position-sensitive detector. The surface neutron source was used to replace the accelerator neutron 
source and moderator to simulate the neutron beams emitted from an ideal neutron aperture. The 
converter screen material was °LiF/ZnS and its specific parameters are listed in Table 2 P", The 
photon-emission spectra of the converter screen and other scintillation parameters were set 
according to previous studies P?! The position-sensitive detector was affixed to the back of the 
converter screen. The single pixel sampling area of the detector is 0.2 x 0.2 mm?, which is 
consistent with the projection of CCD camera pixels on the converter screen. This detector is used 
to record information about the number and position of photons produced from the converter 


3] which generates grayscale images for comparison with the experimental results. Based 


screen 
on this model, the effect of neutron beam characteristics on resolution can be analyzed using an 


arbitrary L/D parameter. 
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Neutron beam 
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Fig.6 Model of the ideal TNR system used in the GEANT4 simulation program. 


The edge-spread function (ESF) method was used to extract the resolutions at different 
locations in the simulated image 41. In this method, an opaque sample with a straight edge is used 
as the image object, and the curve of the grey-level change in the direction of the vertical image 
edge in the image plane is called the ESF curve. The spatial resolution was defined as the distance 
between 10% and 90% of the ESF curve P5. A schematic of the resolution measurement 
simulation model using the ESF method is shown in Fig. 7. The red box in the figure represents 
the imaging sampling area, and the center of the red box is located on the edge of the sample. The 
relative distance between the edge of the sample and the imaging center was defined as de. The 
sample used in the experiment was a | cm-thick rectangular iron block. Different images were 
obtained in the simulations by adjusting the edge positions of the samples. The curves of the grey 
values of the images were obtained from the simulation images, and the resolution of each 
position was calculated using the ESF method. 


Table 2 Optical properties of the converter screen 7! 


Properties Specification 
Material SLiF/ZnS 1:2 
Dopant Ag 
Density 2.6 g/cm? 

Optical photon yield/neutron 160000 
Decay time 1 ps 
Converter screen thickness 400 um 
fa a ea OO 


i 

} 
E 

Distance from the image position to the center 

Converter screen 


Fig.7 Schematic of the resolution of each position measured using the ESF method. 
The simulation results were compared with experimental images to verify the accuracy of the 
simulation. The yellow boxes shown in Fig. 4 (L/D=7.4) are selected as the sampling areas with 
d=1.0 cm and d=5.0 cm. The grey images extracted from the experimental image are presented 
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together with the images obtained from the simulations with the same parameters of de in Fig. 8a. 
The gray curves extracted from these images are compared in Fig. 8b. The results show that the 
simulated curves agree well with the experimental curves. The calculated resolution of the 
simulated image deviated from the experimental data by less than 4.5%. 


(a) Experiment: (b) 1.0 


Simulation result: 
1 cm 


$ cm 
Experiment result: 
aoaia Lom 


. . =-=: Scm 
Simulation: 


ari 5 10 T5 20 25 30 35 40 we 50 
Pixel 
Fig.8 Comparison of experimental and simulation results: (a) grey image of sample edge and (b) corresponding grey 
curve. 

Based on the validated simulation analysis method, the imaging resolutions with L/D=7.4 at 
different positions were first calculated. The simulation results are shown in Fig. 9a. The black 
line represents the spatial resolutions at different positions, as shown in Fig. 7, and the red line 
represents the deviations between the simulated and theoretical resolutions. In Fig. 9a, the 
resolution curve exhibits a symmetrical structure, and the deviations increase linearly with de. To 
further study the influence of de on the resolution under different L/D conditions, images with five 
different L/D ratios were simulated, and resolutions at different imaging positions were obtained. 
Fig. 9b presents the variations in the deviations between the simulation and theoretical calculation 
results at different de positions. The results show an incremental deviation as de increases, which 
gradually narrows as L/D increases. Another interesting result (Fig. 9b) is that the simulation 
results are consistent with the resolutions of the theoretical model at the center of the image (d-=0), 
and the maximum deviation is less than 2.7%. 

The experimental data analysis provided a result similar to that of the Monte Carlo simulation. 
Fig. 9c presents the ESF curves derived from the experimental data for L/D=7.4 and L/D=12. 
Under different L/D conditions, there was a significant deviation between the imaging locations 
d.=1 cm and d.=5 cm. Moreover, the deviation with a smaller L/D ratio is larger. The experimental 
data intuitively show the impact of geometric distortion on the imaging resolution. 
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Fig.9 Resolution analysis of geometric distortion by LSF method: (a) simulation results of relationship between 
resolution and d. with L/D=7.4, (b) variation of deviations between simulation results and theoretical 
calculation results with different L/D ratios, and (c) experimental ESF curve comparison of line-pair edge 
images with different L/D ratios. 
Based on the above analysis, it can be concluded that geometric distortion is an important 
factor that leads to resolution deterioration. Under low-collimation conditions, the resolution 


deterioration caused by geometric distortion increases linearly with de. 

To explore the mechanism underlying this phenomenon, the neutron beam emission angle 
distribution was calculated using a previously built Monte Carlo simulation model. Neutron 
emission angles were recorded on the front surface of the converter screen at different de intervals. 
Fig. 10a presents the calculation results with the de value intervals of 0—1 cm, 1-2 cm, 3—4 cm, 
5—6 cm, 7-8 cm, and 9-10 cm under the condition of L/D=10; Fig. 10b presents the calculation 
results of the de value range of 4.5-5.5 with 5 different L/D ratios. 
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Fig.10 Results of neutron angle distribution with different (a) de values and (b) L/D values. 

The results in Fig. 10 show that when only de increased, the Gaussian widths of the neutron 
emission angle did not change significantly, whereas the average neutron emission angle gradually 
increased. When de is fixed, the average values of neutron emission angle decrease with L/D. 
These results indicate that the degradation of the image resolution is related to the average value 
of the neutron emission angle, which can be derived from the parameter de. 

Based on the above considerations, a theoretical derivation was performed in this study. Fig. 
11 shows the geometric layout of the collimator, object, and converter screen. The offset distance 
dg between the original and actual positions can be calculated based on a simple trigonometric 
relation: 


d = t 6) 


where de is the distance from the object to the imaging center, L is the distance from the aperture 
to the object, and d is the distance from the object to the screen. 
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Fig.11 Schematic of geometrical distortion in neutron radiography. 
An oblique-incidence neutron beam also results in geometric distortion on the converter 
screen. A schematic of the NIL process on the converter screen is shown in Fig. 12. The offset 
distance can be calculated using Eq. 4 as follows: 
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d af oe” (4) 


where d; is the offset of the light output on the converter; ¢ is the thickness of the converter screen 
as the effective photo yield is generated mainly by neutrons that react at the bottom of the screen; 
and L is the distance from the collimator aperture to the converter screen. The offset is 


proportional to ¢ and d and inversely proportional to the L/D ratio. Fig. 13 shows the 


Cc 


simulation results of the light spot generated by a neutron beam with different d , under the 


L/D=10 condition. A phenomenon in which ds moves proportionally with d. can be observed. 
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Fig.13 Light spots generated by the neutron beam vary with de =0 (a), de =5 (b), and de =10 (c). 

Both offsets introduced before blurred the image of the object edge and deteriorated the 
resolution. Geometric distortion is caused by the physical characteristics of the beam emission 
angle. This can also be observed in the traditional TNR and other optical imaging systems °®. In 
contrast to the CTNR system, L/D is usually very large in these imaging areas, and the effect of 
geometric distortion is usually insignificant. 


3.1.2 Effect on converter screen characteristics 


The characteristic parameters of a converter screen include the photon yield and optical 
diffusion response 771. A series of calculations was performed to study the characteristics of the 
converter screen. In the neutron radiography system, the optical diffusion response (6) of the 
converter screen can be calculated as follows !!41: 


h, (%y) =exp[-1X——“] © 


where h 7 (x,y) is the light output distribution induced by a single neutron beam from the 
converter screen. Fig. 14a shows the Monte Carlo simulation results of the photon output behind 
the converter screen, which conform to the Gaussian distribution of Eq. 5. The Gaussian widths of 
the photon distribution with a large L/D ratio were obviously smaller than those with a small L/D 
ratio, whereas the neutron light yield in the spot center decreased with the deterioration of L/D. 
Fig. 14b compares 6 fitting from Fig. 14a and the total neutron light yield together. 5 of the 


converter screen shows a negative relationship with the L/D ratio, which is consistent with the 
analysis results presented in Fig. 14a. However, under different L/D conditions, the total photon 
yield of the neutrons was approximately the same. This indicates that the L/D ratio did not affect 
the light output of the converter screen. Fig. 14c shows the variation curve of the optical diffusion 
response versus screen thickness under different L/D conditions. When L/D=100, the curve shows 
few correlations; however, the other two curves exhibit evident positive correlations. In addition, 
the effect of conversion screen thickness on 6 increases significantly with the decrease in 
collimation ratio. 

In an ideal neutron radiography model, different components are considered independent of 
each other. However, from the analyses above, we can find that under a low L/D condition, a 
change in L/D ratio will affect 5 of the converter screen. This is one of the factors potentially 
leading to the deterioration of the resolution observed in the experiments. When the thickness of 
the converter screen is set to 100 um as used in the experiment, ò is 6.32 um with L/D=10 and 
3.49 um with L/D=100. In this situation, the MTF change caused by 6 is only 0.89% on the 
limiting resolution of devices (1 Ip/mm). This result shows that the change of ò under different 
L/D conditions can affect the image resolution. However, the effect of 5 is not obvious. 
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Fig.14 Different characteristics of converter screen: (a) photon output distribution with different L/D ratios, (b) 
comparison between photon yield and 6 with different L/D ratios, and (c) relationship between 6 and the 


thickness of converter screen. 


3.2 Effect of high-energy neutrons on CTNR system 


3.2.1 Effect on converter screen 


In the CTNR system, a moderator is used to moderate fast neutrons from the accelerator. To 
study the relationship between the neutron energy and spatial resolution, a CTNR system model is 
proposed, as shown in Fig. 15a. An 8 cm-thick lead layer was used as the first layer of the 
moderator and a 5 cm-thick polyethylene layer was used as the second layer to obtain thermal 
neutrons. The reflective layer material was 15 cm-thick graphite. This design can significantly 
improve the thermal neutron flux 8179], The neutron source with an energy of 14 MeV was placed 
on the central axis of the moderator before the moderator, and the beam spot was 1 cm in diameter. 
This model can be used to study the thermal neutron flux, neutron beam energy spectrum, and 
photon yield of neutrons with different energies after fast neutrons pass through different 
moderation materials. Fig. 15b presents the energy spectrum of the neutron beam exiting the 
moderator. The proportion of thermal neutrons in the neutron beam was only 9.13%. After 
moderation, a large number of nonthermal neutrons still existed in the neutron beam. Therefore, 


this model is appropriate for studying the effects of nonthermal neutrons on imaging. 
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Fig.15 Simulation of the CTNR system: (a) GEANT4 model and (b) neutron energy spectrum on the converter 
screen. 

Nonthermal neutrons work with the converter screen through different cross sections. 
Therefore, the photon yield with neutrons of different energies was first calculated (Fig. 16a). The 
photon yield was roughly consistent with the microscopic cross section of °Li. Fig. 16b shows the 
simulation results of the photon yield produced by the neutron beam exiting the moderator. The 
proportion of photons generated by thermal neutrons was 95.79%, whereas that of photons 
generated by nonthermal neutrons was 4.21%. Based on the above results, the effect of 
nonthermal neutrons on the imaging resolution can be neglected in CTNR systems, even with a 
low proportion of thermal neutron beams. 
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Fig.16 Simulation results of photon yield: (a) photon yield of different energy neutrons and (b) light yield distribution 


of neutrons exiting the moderator. 


3.2.2 Effect on CCD 


The shielding of a CTNR system is limited by its geometric dimensions, and it is difficult to 
completely shield high-energy neutrons and their secondary gamma rays. When radiation particles 
reach the CCD chip, they deposit energy that generates electronic noises B’. The noisy images 
obtained under different experimental conditions are shown in Fig. 17. As shown in the figure, the 
effect of white spot noise cannot be neglected, even with shielding. The proportions of noise in the 
neutron radiography images with and without shielding were 2.07% and 7.09%, respectively. 

The presence of radiation noise in images is harmful to the quantitative analysis of the 
resolution. The ESF and MTF methods are based on the numerical analysis of the gray curve in 
the sampling area. Because of the limited sampling area, a few random white spots significantly 
interfered with the grayscale curve, which in turn affected the resolution analysis results. Median 
filtering is sometimes used to eliminate noise in neutron radiography images. However, median 


filtering also leads to an overall decrease in imaging resolution, and the algorithm cannot 


eliminate white spots completely P. The theoretical resolution model does not consider the effect 
of electronic noise and background radiation, which is negligible compared with the radiation 
noise in a short imaging time. A quantitative study of the effect of radiation noise on image 


resolution is introduced in the next section. 


Fig.17 Noisy images under different experimental conditions: (a) before the experiment, (b) without shielding, and (c) 


with shielding. 


4 Optimization of the theoretical resolution model for CTNR 


system 


4.1 Correction of geometric distortion 


The geometric distortion effect deviates from the projection position of the object, which 
leads to the horizontal stretching of the ESF curve. According to the above analysis, the total 
stretch length of the ESF curve is the sum of the spot position offsets Dg and Ds in Eqs. 3 and 4. 
The ESF after geometric distortion can be calculated using Eq. 6: 


(t +d ) 


d 
ESF p ZES F 9+0.8x—*——— 6) 


where ESF, denotes the actual resolution. ESFo is the resolution at d.=0. 0.8, and is the influence 
coefficient defined by the ESF method (the distance between the 10% and 90% grey values of the 
ESF curve). Table 3 presents a comparison of the resolutions of the experimental, simulation, and 
theoretical results. These three results are in good agreement (deviation within 3%). 


Table 3 Comparison of resolutions obtained by different methods 


á Experiment Simulation Calculation Deviation 
(cm) 
1 1.70 1.74 1.70 2.30% 
3 - 2.14 2.08 2.88% 
5 2.39 2.50 2.45 1.97% 
7 - 2.84 2.83 0.35% 


In the measurements, the MTF corresponding to a line pair was calculated using Eq. 1, and 
the input MTF (Mn) was fixed. Owing to geometric distortion, the gradient of the curve 
corresponding to the grey gradient of the line pair changed. However, this did not affect inputs Imax 
and Imin in the measurement of each line pair. The variety of MTF caused by the geometric 
distortion can be expressed as 

MT F (Inas nin) 


— 
Mr Foo ~ Emax Fri) 


where MTF» is the theoretical MTF, which is not affected by geometric distortion, and 
MT F b is the MTF measured experimentally. The image of a slit was obtained by 
subtracting the ESF curves of the two edges of the slit. Because of the narrow width of the slit, the 


two curves are assumed to be identical. I -I is proportional to the 


ma x mi n 
slope of the ESF curve. The slope is inversely proportional to the ESF resolution of the curve. 
According to Eq. 7, the measured MT F p can be expressed as follows: 


MT F ESF 
——— m == R) 
MT F o ES F p 


where ESF, denotes the ESF resolution of the line pair. ESFo is the ESF resolution at d=0. The 
relationship between ESF, and ESF» is obtained using Eq. 6. When the corresponding resolution 
of ESFo and the theoretical resolution are assumed to be the same, the MTF at the line pair can be 
calculated as: 


M Fo 
MT F , = O 
(140.8 — <) 


This equation was used to analyze the experimental data for calculating the MTF considering 
geometric distortion. The results were compared with those without correction and are shown in 
Table 4. After considering the influence of the geometric distortion, the maximum deviation was 
reduced from 23.1% to 8.82%. This indicates that geometric distortion is an important factor 
affecting the resolution of the CTNR system. 


Table 4 Comparison of resolutions between experiment and the correction model considering geometric distortion 


Resolution of Resolution of Deviation between Resolution of correction Deviation between 
L/D traditional model ; experiment and model considering experiment and correction 
experiment /mm a eee 
/mm traditional model geometric distortion /mm model 
6.6 1.77 89 6.35% 1.81 441% 
7.4 1.51 .73 12.72% 1.61 7.45% 
8.8 1.27 .48 14.19% 1.36 8.82% 
10 1.11 .28 13.28% 1.23 4.06% 
12 0.91 12 18.75% 1.03 8.04% 


4.2 Correction of the radiation noise 


The effect of radiation noise can be regarded as an independent factor, considering the 
irregularity in the scattering of fast neutrons. When the effect of noise is considered, the total MTF 
can be calculated as: 


MT F ToT (u )= 
. d 2 ; As 
sinc u -exp(-—a (6u sinc u -MT F 
Gie pepr (6u 4) <e 
(10) 
where MT F ToT is the total MTF value containing noise and MTF, is the MTF of 
the noise. 


This was because the noise distribution was random. Eq. | can be used to derive MTF;. 
Within a fixed imaging area, the number of samplings included in the area is n, the proportion of 
noise in the image is o, and the average gray value of noise is g. When noise is randomly added to 
the image, the increase in the gray value of each column is ong. Hence, Laax = Imax + 
ong and T 


min 


= Imn tong . The MTF, calculation formula is expressed as: 


M out  _ (Gnatong)-Cnin* 0M8)1/ [Cmax 808)+Tnin tOn] 1 an 
M | (nax-Imin)/ Anax in) 1+20ng/ Anax Inin) 


in 


MTF,= 


The grey value of the noise in an image is constant in a certain neutron radiography system. For 
single imaging, Imax and Imin are regarded as constants. Therefore, ng/(Lnmax tmin) in Eq. 5 can be 
replaced by the constant r. Eq. 11 can be simplified as: 

MTF, =— 


1+20r 


(12) 

Eq. 12 shows that the proportion of noise directly affects the modulation of the image. The MTF 
was negatively correlated with the proportion of noise. Experimental verification of Eq. 12 has 
been performed in detail in the literature B”. 

In this study, white spot noise parameters were extracted from the experimental images to 
correct the theoretical calculation results. First, the white spot noise was recognized according to 
the abnormal gray gradient from the experimental image in Fig. 4. The average gray noise values 
and proportion of noise were calculated for each image. Subsequently, MTF, for each image was 
calculated using Eq. 12. Table 5 lists the corrected theoretical and experimental results. After 
considering both the effects of geometric distortion and radiation noise, the maximum deviation 
between the theoretical model and the original experimental results was reduced to 4.22%. A small 
systematic deviation between the revised theoretical model and the experimental values was still 
observed. This is mainly due to the interference of non-fixed device factors, such as neutron 


scattering by the experimental sample. 
Table 5 Comparison of resolutions between experiment and corrected model considering geometric 
distortion and radiation noise 


: Resolution of correction Deviation 
Resolution of ; na 
L/D “ae Resolution of model considering between 
traditional : MTF, mee i , 
experiment /mm geometric distortion and experiment and 
model /mm bie : Í 
radiation noise /mm correction model 
6.6 1.77 1.89 0.851 1.86 1.07% 
7.4 1.51 1.73 0.805 1.67 3.59% 
8.8 1.27 1.48 0.813 1.42 4.22% 
10 1.11 1.28 0.909 1.24 3.13% 
12 0.91 1.12 0.791 1.09 2.75% 


5 Conclusion 


A significant deviation in the resolution between the experiment and the traditional theoretical 
resolution model was observed in a neutron radiography experiment performed by the Hefei 
Institute of Physical Science. To explain this phenomenon, the effects of a low L/D ratio and 
high-energy neutrons on imaging resolution were analyzed using Monte Carlo simulations. The 
results indicated that the deviation was primarily caused by geometric distortion and radiation 
noise. Additionally, through theoretical model optimization, the resolution deviation was 
effectively reduced by considering the imaging position and introducing a noise factor. The 
following conclusions were drawn. 

1) When the imaging position is far from the imaging center, geometric distortion leads to the 
degradation of image resolution in the CTNR system. The geometric distortion becomes 
significant with a decrease in L/D. 

2) The optical diffusion response of the converter screen was not independent of L/D under 
low L/D conditions. However, a change in 6 has little effect on the resolution in a thin converter 
screen. 

3) The effect of nonthermal neutrons on the spatial resolution can be neglected in the CTNR 
system. However, the radiation noise induced by high-energy neutrons and secondary y -rays will 


lead to a resolution degradation. 

4) Considering the effects of geometric distortion and radiation noise, the modified 
theoretical model agreed well with the experimental results. The maximum deviation decreased 
from 23.1% to 4.22%. 

This study analyzed the factors affecting imaging resolution in a low-L/D environment, 
which is the main difference between the CTNR system and a traditional TNR system. Currently, 
resolution degradation caused by geometric distortion and radiation noise has not been 
fundamentally solved. The design improvement of the moderating collimator and performance 
optimization of imaging systems are key directions for future research. 
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